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Abstract 

The origin of the most energetic cosmic ray events, with E > 10 20 eV, remains myste- 
rious. One possibility is that they are produced in the decay of very massive, long-lived 
particles. It has been suggested that these so-called "top-down scenarios" can be tested 
by searching for ultrarelativistic neutralinos, which would be produced copiously if super- 
particles exist at or near the TeV scale. In this paper we present a detailed analysis of the 
interactions of such neutralinos with ordinary matter. To this end we compute several new 
contributions to the total interaction cross section; in particular, the case of higgsino-like 
neutralinos is treated for the first time. We also carefully solve the transport equations. 
We show that a semi-analytical solution that has been used in the literature to treat the 
somewhat analogous propagation of r neutrinos leads to large errors; we therefore use a 
straightforward numerical method to solve these integro-differential equations. 



1 Introduction 



Experiments like Fly's Eye [JJ and AGASA jZ] have shown the existence of cosmic rays with 
energy E > 10 20 eV, the so-called ultra high energy (UHE) component of the cosmic ray 
spectrum. This raises two questions: how were these UHE particles produced, and how did 
they manage to reach Earth? 

In order to accelerate ultra-relativistic charged particles to energy E, one needs a magnetic 
field with strength B extending at least over the confinement radius L oc Z/(B ■ E), where Z 
is the charge of the particle [3]. The problem is that very few, if any, objects in the Universe 
can be said with confidence to have a sufficiently large B ■ L to accelerate protons to 10 20 eV, 
once synchrotron losses in high B— fields have been taken into account. 

The second problem is that the arrival directions of UHE events seem to be distributed 
more or less homogeneously. Since protons of this energy would not be deflected much by the 
magnetic fields in our galaxy, this excludes one or a few local point sources. A large group 
of very distant sources, e.g. related to active galactic nuclei, would indeed yield an essentially 
homogeneous distribution. Unfortunately particles produced at such very large distances should 
not be able to reach us without losing much of their energy. In particular, as pointed out in 
refs. jlj, protons with E > 5 • 10 19 eV would lose their energy through inelastic scattering on the 
photons of the cosmic microwave background (CMB); this is known as GZK effect or GZK cut- 
off. The same target also depletes the flux of photons with E > 10 15 eV, the most important 
reaction here being the production of e + e~ pairs. Finally, heavier ions would suffer break- 
up reactions on CMB photons. The upshot is that no known UHE particle which interacts 
high in the atmosphere (as required by observation) would be able to travel over distances 
exceeding ~ 50 Mpc [3]. Intergalactic magnetic fields should not be able to randomize the 
arrival directions of UHE particles over such distances, i.e. they should still (more or less) point 
back to their sources. However, there are no known sources within this radius in the directions 
of these UHE events. 

This has lead to various exotic proposals [HIE!- In particular, it has been suggested [HI Ej 
that UHE events originate from the decay of very massive, long-lived particles X with mass 
Mx ^ 10 12 GeV. These could be particles associated with a Grand Unified theory, which are 
stabilized by being bound in topological defects [HI EJ, or free particles with extremely small 
couplings to ordinary matter [$,. In either case, they would have been produced in the very 
early universe, just after the end of inflation [SI El- The GZK effect could then be circumvented, 
if most of these decays occur within one GZK interaction length, for example in the halo of our 
own galaxy. The decay of an X particle triggers a parton cascade, followed by hadronization 
and the decay of unstable particles [TO*! ITT]. 

One potential problem of this scenario is that it tends to predict a higher flux for UHE 
photons than protons, at least at the point of X decay. For E <C Mx this is due to the fact 
that fragmentation produces more pions than protons, and neutral pions decay into photons; 
for E ~ the direct emission of photons in the early stages of the parton cascade also 

plays a role. This is problematic, since experiments indicate that most UHE events are caused 
by protons (or heavier nuclei) rather than photons [12]. However, the statistics at post-GZK 
energies is still quite poor, and our knowledge of the details of interactions at E > 10 19 eV 
(corresponding to ^fs > 10 5 GeV) still leaves much to be desired. Propagation effects might 
also modify the proton to photon flux ratio. Finally, protons from conventional ( "bottom-up" ) 
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sources might contribute significantly (although not dominantly) to the observed UHE events, 
thereby increasing the proton to photon ratio ^3j- Given that other explanations also have 
their problems top-down models should still be considered viable. 

These models predict an UHE neutrino flux that is even higher than the photon flux, since 
fragmentation produces more charged than neutral pions, and more W bosons than photons are 
emitted during the early stage of the parton cascade. These models can therefore also be tested 
[14] by looking for neutrinos with E > 10 5 GeV, where the atmospheric neutrino background 
becomes negligible. However, bottom-up models also predict a substantial UHE neutrino flux, 
e.g. due to the GZK effect itself. Detailed analyses of the spectra of neutrino events would 
therefore be required to distinguish between top-down and bottom-up models. 

There is, however, a potential "smoking gun" signature for top-down models. The hierarchy 
between the scale Mx ^ 10 12 GeV and the weak scale will in general only be stable against 
radiative corrections in the presence of weak-scale supersymmetry ^H] ■ If R parity is conserved, 
the lightest superparticle (LSP) is stable, and will also be produced copiously in the decay of 
X particles [To^ lllj. In most supersymmetric models, the best LSP candidate (in the visible 
sector) is the lightest neutralino Xi [13 • Early estimates [TH] indicated that the cross section 
for neutralino interactions with ordinary matter is significantly smaller than that of neutrinos. 
This lead to the suggestion [TZJ [TBI EH] to search for UHE neutralinos by using the Earth as a 
filter: for some range of energies the neutrino spectrum would get depleted in the Earth, but 
neutralinos would still come through, and might be detectable by future experiments. 

Clearly a good understanding of the interactions of UHE neutralinos with matter is manda- 
tory in order to analyze the viability of this signal. This is the topic of our paper. The 
propagation of neutralinos through the Earth was treated using a simple Monte Carlo model 
in ref.^7], whereas ref.^Hl ignored all neutralinos that had any interactions in the Earth (as 
opposed to in the detection medium). This is "overkill", since for models with conserved R 
parity, any neutralino interaction will eventually result in another LSP, i.e. another neutralino, 
albeit with reduced energy. In this regard neutralinos are like r neutrinos, which are also 
regenerated after interacting with ordinary matter |2~T]. In order to treat this effect, one 
has to know the interaction cross section differential in the scaling variable y = 1 — E ont /E[ n , 
where E[ n and E out are the energies of the incoming and outgoing neutralino, respectively. The 
calculation of this cross section is described in Sec. 2, including for the first time the case of 
higgsino-like neutralino. In Sec. 3 we describe neutralino propagation through matter by means 
of transport equations, similar to the ones used to describe v T propagation j2I] . We find that 
an iterative solution of this equation along the lines of refs. |22| l2*T| l2"3~] is not applicable in our 
case, since it leads to a significant violation of the conservation of the total neutralino flux; this 
flux conservation is a direct consequence of the regeneration mechanism described above. We 
therefore use a straightforward numerical integration of the transport equations. Finally, Sec. 4 
contains a brief summary and some conclusions. The calculation of a new, usually subdominant 
contribution to the neutralino nucleon scattering cross section is described in an Appendix. 

2 Calculation of cross sections 

This Section deals with the calculation of the total and differential cross section for the scat- 
tering of ultra-relativistic neutralinos Xi on nuclei. This process dominates by far the to- 
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Figure 1: s— channel Feynman diagram for XiQi scattering, where g~j is a virtual squark and the 
symbol X stands for all allowed final states. The arrows indicate the flow of baryon number. 
An analogous diagram exists for XiQi scattering, with reversed arrows. 



tal interaction of neutralinos with matter, except for a very narrow range of energies around 
E^o = (mf — m 2 - )/(2m e ), where on-shell selectrons e can be produced in x°e scattering [19|. In 
the following two subsections we will discuss s— channel and t— channel scattering, respectively; 
we will see that the former (latter) process dominates for bino- (higgsino-)like LSP. 

2.1 s— channel contribution 

The s— channel contribution to neutralino-quark scattering, XiQi ~^ X, is described by the 
Feynman diagram shown in Fig. 1. The total partonic cross section can be written as follows: 

I 1 I v Hi' Hi Hi 

where s is the partonic center-of-mass (cm.) energy, is the squark mass, P* = (s — 
m 2 ~ )/{2ys) is the cm. 3-momentum of the incoming particles, is the total decay width of 
the squark and r(g~j — > qt + xi) is the partial % — > + Xi decay width. The same expression 
also holds for XiQi scattering. 

The s— channel contribution to the total Xi~ nucleon scattering cross section can be eval- 
uated from Eq. flj) by convoluting with the appropriate (anti-)quark distribution function and 
summing over flavors. We use the narrow width approximation, 
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(s - m 2 ) + m 2 r 2 rnT 
The convolutions then collapse to simple products [T%] : 



<5(s - m 2 ) . (2) 
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x= * * (4) 



2m nE 
3 



m 



where E- in is the energy of the incident neutralino in the rest frame of the nucleon, and = 
(nip + m n )/2 is the nucleon mass. In our numerical evaluation we use Q 2 = m| as momentum 
scale in the (anti-) quark distribution functions. For simplicity we assume equal masses m q for 
the L and R squarks of a given flavor. In this case left- and right-handed couplings contribute 
symmetrically, as shown in Eq.©. 1 Our simplified treatment is sufficient as long as squark 
masses remain free parameters; note also that most SUSY models predict small mass splittings 
between squarks, at least for the first two generations [TH] . Finally, the couplings appearing in 
Eq.(jnj) are given by: 

a qL = V2g 2 \T3, q N 12 + tan ^ jy u ^ f or q = Uy d,s,c; 
a qR = \Plg 2 tan OwQqNu for q = u,d,s,c; 

9 2<m b M f ,s 
a b L ,R - a d L ,R + —7E 7 iV 13 • yO) 

V2rriw cos p 



Here, Ny are the entries of the neutralino mixing matrix in the notation of ref . |24j . g 2 is the 
SU{2) coupling constant, Ow is the weak mixing angle, T 3u = —T 3 d = 1/2 is the weak isospin, 
Q q is the electric charge of quark q in units of the proton charge, mw is the mass of the W 
boson, tan (3 is the ratio of the Higgs vacuum expectation values, and mb is the mass of the 
bottom quark. We ignore the masses of quarks of the first and second generation in these 
couplings, i.e. we use identical couplings for up and charm quarks, as well as for down and 
strange quarks. We do include contributions oc m;, to the couplings of the bottom quark [2"4] . 

Note that we do not include the contribution from top (s)quarks in Eq.(j3J). Since the top 
quark may not be much lighter than their superpartners, it is more appropriate to treat t 
production through the 2-^2 scattering reactions 

X° l9 -> it, tt . (6) 

We evaluated the corresponding cross section, but found it to be subdominant in all scenarios 
we considered; this is not very surprising, since it is of higher order in the strong coupling than 
the cross section (J3J). We therefore delegate the discussion of the reaction (JBJ) to the Appendix. 2 
In order to define potentially realistic supersymmetric scenarios, we work in the framework 
of minimal supergravity (mSUGRA) fTS], which is the simplest supersymmetric model generi- 
cally yielding a stable neutralino as LSP (assuming that the gravitino is heavier than the Xi)- 
Here one postulates universal sfermion (mo) and gaugino masses (mi/2) at the scale of Grand 
Unification. We use the public code Softsusy [22] to calculate spectra for three representative 
scenarios yielding a bino-like state as LSP, as defined in Table ^ For simplicity we fix the 
trilinear interaction parameter A = and tan/? = 2. We choose rather small or moderate 
values of the scalar and gaugino mass parameters; larger masses yield smaller cross sections, 
and hence smaller effects from propagation through Earth. 3 

general a qL only contribute to (}l exchange, while a qR contribute to exchange. 
2 One should not add the cross section from x?<7 - * QQ to the cross section ©, as is done in ref. [lfi|. 
The reason is that the result © already includes contributions from gx_i scattering in the leading logarithmic 
approximation, via g — > qq splitting, which contributes to the Q 2 dependence of the quark distribution functions. 
Simply adding both contributions therefore leads to double counting. Instead, gxi scattering should be treated 
as part of higher-order corrections to inclusive q production described by Eq. @ , provided the quark q can be 
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mSUGRA scenarios 


Scenario 


m 


mi/2 


Xi 


m d L 


Dl 


80 


150 


63 


365 


D2 


150 


250 


104 


582 


D3 


250 


450 


189 


992 



Table 1: mSUGRA scenarios for tan/3 = 2, higgsino mass parameter fi < and A = 0. m 
and mi/2 are the universal scalar and gaugino mass parameters, respectively. m^ L is the mass 
of the SU(2) doublet down squark, but all other squark masses have quite similar values. All 
masses are in GeV. 



Results for the corresponding total cross sections for scattering on nucleons, with and with- 
out the contribution from bottom (s)quarks, are presented in Fig. [21 Here and in the following 
figures we use the CTEQ6 parameterization [27 j of the parton distribution functions, averaged 
over proton and neutron targets; in the absence of large mass splitting between u nd d squarks, 
the cross sections for scattering on protons and neutrons are very similar. We see that the 
contribution of bottom (s)quarks only reaches 10% even for large squark masses and large Xi 
energy. In this case one is probing the b— quark distribution function at small Bjorken— x, see 
Eq.pjl. and rather high Q 2 , where it is comparable to the distribution functions of light quarks. 
However, the small hypercharge of bottom quarks suppresses their couplings to bino-like neu- 
tralinos, relative to those of up and charm quarks; see Eq.(j3J). 

Over most of the parameter space of mSUGRA the lightest neutralino is bino dominated 
[l3|l26j: this is also true for the three scenarios of Tabled Since the bino coupling to a sfermion 
is proportional to the hypercharge of that sfermion, a bino-like Xi couples predominantly to 
SU(2) singlet, "right-handed" squarks, whose hypercharges are two (for cLr) or four (for ur) 
times larger than that of the SU(2) doublet squarks; see Eqs.(jHJ). The total s— channel cross 
section is therefore dominated by the production of 577(2) singlet squarks. Since these squarks 
do not couple to SU{2) gauginos, they will decay directly into q + Xi, if the gluino is heavier 
than these squarks [2E]; i n mSUGRA this corresponds to m < mi/2- In this case we can 
approximate the total s— channel contribution as production of on-shell squarks which decay 
back into q + Xi fi na l states. 

This greatly simplifies the calculation of the cross section differential in the scaling variable 
y = 1 — E out /E in , where E in and E out are the incoming and outgoing Xi energy in the nucleon 
rest frame. The crucial observation is that squark decays are isotropic in the squark rest frame, 
which implies 



dr. a 



tot 



(7) 

dcosfl* 2 ' y ! 

where 8* is the angle between the ingoing and outgoing Xi i n this frame. In order to get the 

assigned a partonic distribution function in the nucleon. 

3 Within mSUGRA, the combination of small tan/3 and small soft breaking masses is excluded by Higgs 
searches at LEP. Note, however, that tan/?, or indeed the structure of the Higgs sector as a whole, plays little 
role in our analysis. Scenarios with rather light squarks and gauginos are still allowed in other supersymmetric 
models. The squark masses in scenarios D2 and D3 are well above the lower bound derived in |2(ij in the 
framework of mSUGRA. 
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Figure 2: The total cross section of the s— channel contributions to — nucleon scattering for 
four (dotted curves) and five (solid curves) active quark flavors. The super symmetric scenarios 
Dl, D2 and D3 are defined in Tabled 



final expression for the y distribution, we have to boost from the cm. system into the rest 
frame of the nucleon. This yields a flat distribution, 



cb s al ot 



where 



(8) 



1/max 1 



m? 

Xi . 

2 ' 



Vmin = 0. (9) 

In the first Eq.@ we have used s = m~ for on-shell squark production. Forward scattering 
in the squark rest frame leads to E out = E- m , independent of the details of the kinematics; the 
lower limit for y is therefore always zero. On the other hand, y = 1, which requires E out = 0, is 
possible only for m^o — > 0. Since in mSUGRA m? m§, |/ max is indeed quite close to unity. 
An LSP will then lose on average about half its energy if it undergoes s— channel scattering on 
a nucleon. 



2.2 t— channel contribution 

We now turn to the calculation of the t— channel contribution to the LSP nucleon scattering 
cross section. As shown in Fig. El there are both W (charged current) and Z exchange (neutral 
current) diagrams; the former produce an outgoing chargino, while the latter have one of the 
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Figure 3: Feynman diagrams for t— channel neutralino-nucleon scattering. The arrows indicate 
the flow of baryon number, and p± ; 2, ^1,2 and q refer to the four-momenta of the initial, final 
and exchanged particles, respectively. 

four neutralinos in the final state. There are additional diagrams with antiquarks in both 
initial and final state, as well as charged current diagrams producing a negative chargino from 
scattering off a d or u quark. Table [2] shows all possible initial and final states, where we again 
include five active flavors of quarks in the nucleon. 

The partonic total cross section can be obtained by integrating over the scattering angle 8* 
in the center-of-mass frame; convolution with the relevant quark distribution functions then 
yields the t— channel contribution to the Xi~ nucleon scattering cross section: 



<io) 



where s = 2xEi n mN+m 2 , and P* and P 2 * are t ne three-momenta of the incoming and outgoing 
particles in the cm. system, respectively. As usually done in deep-inelastic scattering, we 
identify the scale Q 2 in the quark distribution functions with the absolute value of the four- 
momentum t exchanged between the participating partons. This causes a minor difficulty: 
forward scattering (cos#* = +1) leads to Q 2 = t = 0, where the parton distribution functions 
are not defined. Demanding Q 2 > Q m i n therefore leads to a restriction on the phase space 
integration: 

COsCax=l-7 -, (11) 

[s-ml J (s-mU 

where Xout stands for the chargino or neutralino in the final state. This in turn affects the lower 
bound on the momentum fraction x: 



\ (Qmin + ™jL - m|o) + yj \ (<3min + ^ + m| ? ) ~ ^| out ^|o 
Xmin = op ~ ' (12) 

where E- m is again the energy of the incoming neutralino in the nucleon rest frame. The total 
cross section depends only very weakly on the choice of Q min . The W and Z propagators 
become independent of Q 2 = \t\ once Q 2 <C m^. Moreover, as we will see shortly, the total 
cross section is dominated by the production of Xout 7^ Xi-> where Q min m \ mt om y has a very 
small effect on x min . In our numerical calculations we fixed Q min = 1 GeV 2 . 
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electroweak gauge boson in the t— channel, together with the number of different final states. 
The chargino index k runs from 1 to 2 and the neutralino index / from 1 to 4. 



The squared matrix elements \M\ 2 for the charged current reactions are given by 
1 



Em 



A 



spins 



2\C^\ 2 Pl -p 2 h-k 2 + 2\C? k \ z Pl -k 2 p 2 -k 1 



iR |2„ 



for the first four cases of Table |2l and 



si 



spins 



[2 1 Ctk \ 2 Pi ■ k 2 p 2 ■ h + 2 1 C* | ^ -p 2 fci • fc 2 



(13) 



(14) 



for the last six cases of Table 121 note that the results (|13|) and (|14|) differ by the exchange of 
the left- and right-handed W^XiXk couplings. These couplings are given by [21] 



N 12 V* kl - ^n u v; 2 , 



N* 2 U kl + -y=N* 3 U k2 . 



(15) 



Here, the index k runs from 1 to 2, W and V denote the chargino mixing matrices, m^± is the 

Xk 

mass of the outgoing chargino, and t = —Q 2 = (&i — Pi) 2 ', the four-momenta have been defined 
in Fig. El In Eqs.(|13|) and (|T4*|) we have ignored quark flavor mixing, i.e. we replaced the quark 
mixing matrix by the unit matrix. 

Due to the Majorana nature of the neutralinos, neutral current reactions on a quark and 
antiquark are described by the same matrix element: 



spins 



2 gU\gl\ 2 + K 



-,9 I 



-IV 



m 2 w : 



(\N l L 1 \ 2 + \Ng\ 2 ) { Pl -p 2 k 1 -k 2 +p v k 2 p 2 -k 1 ) 



(N&Ng* + N^Nt) p 2 -k 2 m^m^ 



(16) 



with 



c w 



R 
II 

9l 
9 q n 



cos 9w , 

I (~N l3 N* 3 + N l4 N* 14 ) 



■3,9 



Q q sin 2 W . 



Q g sin 2 9 W 



(17) 



The index / in Eqs.([TBjl and (|T7|) runs from 1 to 4, and Q g and T 3 (? have been introduced in 
Eqs.©. 



Scenario 






m x° 3 






m-+ 

X?. 


GF [%} 


D2 


104 


206 


468 


477 


206 


476 


99.4 


D3 


189 


367 


800 


805 


367 


805 


99.8 


HI 


125 


137 


742 


801 


130 


742 


1.2 


H2 


300 


310 


940 


970 


303 


970 


1.1 



Table 3: Neutralino and chargino masses (in GeV) as well as the gaugino fraction GF = 
\N U \ 2 + \N 12 \ 2 for the four scenarios discussed in the text. Scenarios D2 and D3, where Xi 
is bino-like, have already been introduced in Table ^ whereas scenarios HI and H2 describe 
higgsino-like xi', n °t e that the higgsino fraction is 100% — GF. 



As mentioned earlier, all neutralino and chargino states can be produced in t— channel 
Xi~ nucleon scattering. Table El lists the masses of these particles for four scenarios. The 
first two have already been introduced in Table Q and describe bino-dominated xi states, 
whereas in scenarios HI and H2 Xi is dominated by its higgsino components. Eqs. (jl5j) and 
(fT7|) show that bino-like states, which have \Ni 2 \, [ iV" 13 1 , |AT 14 | 1, have suppressed couplings 
to gauge bosons. We therefore expect the t— channel contributions to be subdominant for the 
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scenarios with bino-like Xv This is borne out by Fig. which shows the s— and t— channel 
contributions to the total x\~ nucleon scattering cross section for the three scenarios of Tabled 
The t— channel contribution dominates only at low energies, below the threshold for squark 
production; however, there the cross section is in any case very small, and any possible signal 
from UHE LSPs will be masked by the much larger neutrino signal. At higher energies the 
s— channel contribution is at least a factor of 30 larger than those from all t— channel diagrams. 
Note that the latter also decrease quickly with increasing sparticle mass scale, i.e. when going 
from scenario Dl over D2 to D3. The reason is that increasing the gaugino mass parameter 
mx/2 in mSUGRA not only increases the bino mass, via the condition of electroweak symmetry 
breaking it also increases the absolute value of the higgsino mass. Both effects tend to reduce 
gaugino-higgsino mixing, and hence the couplings of Xi to gauge bosons. These results show 
that one can indeed ignore all t— channel contributions for bino-like LSPs, as done in refs.|17[ 

HE]. 
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Figure 4: The total Xi~ nucleon scattering cross section for the three scenarios of Table ^ 
where Xi is bino-dominated. At large energies the upper (lower) curves show the contribution 
from s— (t— )channel diagrams. Interference between these diagrams is negligible. 



In contrast, Eqs.(j5J) show that a higgsino-like x\i with |iVn|, \N i2 \ <C 1, has very small cou- 
plings to light quarks, which suppresses the s— channel contribution, whereas | iVi3 1 , \N\&\ are 
sizable and lead to unsuppressed couplings to gauge bosons. We therefore expect the t— channel 
contributions to dominate the total x5 — nucleon scattering cross section for higgsino-like Xv 
Fig. shows that this is indeed the case. The solid and dotted curves in this figure show the 
total t— and s— channel contributions in scenarios HI and H2 of Table El where the t— channel 
contributions now correspond to the upper curves. Here we used squark masses of around 1.9 
TeV for scenario H2, and about 1 TeV for HI; note that large scalar masses, rriq 3> are 



10 



required in mSUGRA if the LSP is to be higgsino-like ^H]- For comparison we also show the 
total cross section for neutrino-nucleon scattering [2H] (dot-dashed curve). Since neutrinos are 
also SU(2) doublets, their cross section is similar to that of higgsino-like neutralinos. However, 
some differences remain even at high energies. The reason is that the quark distribution func- 
tions appearing in the expression (jlUj) for the total t— channel contribution to the cross section 
peak strongly at small x. As a result, the partonic cms energy Vs is often not much larger than 
the masses of the relevant charginos and neutralinos, leading to some suppression of the LSP 
cross section relative to the neutrino cross section. We will come back to this point shortly. 




le+05 



le+06 le+07 

E^[GeV] 



le+08 



le+09 



Figure 5: Comparison of total cross sections for neutrinos (dashed-dotted curve) and higgsino- 
and bino-like Xv F° r scenarios HI and H2 (see Table |3J) with higgsino-like LSP, t— and 
s— channel contributions are shown by the upper and lower curves, respectively; for scenarios 
D2 and D3, only the dominant s— channel contributions are shown. 



The dashed and dash-doubledotted curves in Fig. El show again the results for scenarios 
D2 and D3, where we have omitted the subdominant t— channel contributions for clarity. We 
see that these scenarios lead to somewhat smaller total cross sections than the scenarios with 
higgsino-like LSP. The s— channel partonic subprocesses occur at first order in the weak cou- 
pling, whereas the t— channel diagrams are of second order. However, this relative enhancement 
of the s— channel contributions is over-compensated by the fact that the relevant scale in the 
partonic s— channel cross section is the squark mass, while for the t— channel cross sections 
it is given by the mass of the exchanged gauge boson. On the other hand, comparison with 
Fig.0]shows that for scenario Dl with relatively light squarks, the (s— channel dominated) cross 
section is about the same as the (t— channel dominated) cross section for higgsino-like LSP. 

Scenarios with higgsino-like Xi have higgsino mass parameter significantly smaller than 
the SU(2)l and U(l)y gaugino mass parameters |M 2 | and |Mi|. As a result, the second 
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neutralino X2 an d lighter chargino Xi w iU a l so be higgsino-like, but x° and xf w iH be 
gaugino-like. Since gauge bosons can only couple to two higgsinos (or, in case of W ± , to two 
SU{2) L gauginos), the production of these heavier neutralino and chargino states is strongly 
suppressed, as shown in Fig. El This figure shows that the neutral current reaction with x\ m 
the final state is also suppressed. This is due to the fact that |-/V 13 | ~ |AT 14 | for higgsino-like 
LSP, leading to a strong cancellation of the ZxiXi couplings in Eqs. ([T?jl . The total neutral and 
charged current contributions are therefore dominated by X2 an d Xi production, respectively. 
The Xi cross section is about two times larger, since xt an d Xi contribute about equally. 
All three states are essentially SU{2) doublets, and are produced with nearly pure vector 
coupling, leading to nearly identical cross sections for xf production on a quark or antiquark, 
see Eqs.jni), d- 4 




Figure 6: Separate contributions from the four neutralinos and two charginos to the total 
t— channel cross section in scenario H2 of Table El 



We now turn to a discussion of the y— distribution of the t— channel contribution. The 
scaling variable y is related to the scattering angle via 

2 

e — 772 ~ 

y=(l-cos0*) (18) 

The cross section differential in y can therefore be calculated from Eq. (jlO|) by a substitution of 

4 In fact, after summing over quark and antiquark contributions, at high energies the cross sections for \t 
and Xi production are nearly the same even if left- and right-handed couplings of Eas. (|15|) differ from 

each other. The reason is that dominant contributions come from small values of x, where the distribution 
functions for quarks and antiquarks are essentially identical. For the same reason the cross section for the 
scattering of neutrinos and antineutrino becomes identical at high energies |3(J1 129j . 
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variables, and inverting the order of integration: 



Eqs.()18p and (|TT|) imply that for given x, 



\M\ 



16tt(s - m 2 0/ 

Xi 



(19) 



y e 



2 ■ 



2xE- m rriN 



Xout 



(20) 



Using s = m? + 2xE- m mN, this gives for the lower bound on the x— integration in Eq. (fITJ|) : 

Xi 



i 



max 



D 2 m- 

Vmin '""Xout <2 

-m-o 



y 1-2/ 



(21) 



Since Qmin = 1 GeV <^ m , for most values of y the lower bound on x is determined by the 
second term in Eq.(|21|); the restriction on Q 2 determines x_ only for very small values of y. 




Figure 7: Differential neutral and charged current t— channel cross sections da/dy for LSP- 
nucleon scattering in scenario H2 (black curves), and for neutrino-nucleon scattering (red or 
grey curves), for three different energies. The charged current contributions are given by the 
upper curve for a fixed energy. 



Fig. El shows the differential cross section da/dy for scenario H2 and three different energies. 
The charged and neutral current contributions are shown separately. They have very similar 
shapes at all energies. This is not surprising, since in the pure higgsino limit, both contributions 
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are from vector-like couplings with equal strength. Moreover, Table El shows that the masses 
of X2 an d Xi are ver y close to each other; we saw in Fig. El that the production of these two 
particles dominates the total t— channel contribution. 

Fig. [7| also shows the differential cross section for neutrino-nucleon scattering at the same 
energies. While neutrino and neutralino cross sections are quite similar for y <C 1, the latter 
fall much more quickly with increasing y. This is largely due to the lower bound (|2"T|) on the 
integral over x in Eq. (jl9|) . which increases quickly with increasing y. Notice that in the region 
of small y, which dominates the total cross section, x_ is given by the difference between the 
masses of the incoming and outgoing \ particles. Table El shows that this difference is smaller 
for scenario H2, which has much larger masses for the higgsino-like states. This explains why 
the difference between the total cross sections in scenarios HI and H2 depicted in Fig. El is 
smaller than that between the cross section for HI and the neutrino cross section. Scenario H2 
still gives a smaller cross section since the mass terms enter the squared matrix elements p3jl . 
()14|) and (j!6|) with negative sign. The same effect also explains why the ratio of charged and 
neutral current contributions falls slightly with increasing y: the fact that Xi 1S slightly lighter 
than X2 is more important at small y. In case of neutrinos, y — > 1 remains possible at very 
small x. However, increasing y also means increasing the squared four-momentum exchange \t\ 
for fixed s; this explains the decrease of the neutrino cross sections. 

Up to this point we have not considered the interference between the s— and t— channel 
contributions. This is suppressed by two effects. First, it can only occur for outgoing x\ 
state; we saw in Fig. El that this final state does not contribute much to the total t— channel 
cross section. Secondly, the interference between the t— channel gauge boson and s— channel 
squark propagators vanishes at I = m§, where the latter is largest in absolute size, but purely 
imaginary. Interference can thus only occur for off-shell squark exchange, and is hence of higher 
order in perturbation theory. These interference contributions can therefore safely be neglected. 



3 Solution of the transport equation 

We are now ready to discuss the propagation of UHE neutralino LSPs through matter. This 
is described by so-called transport equations, which have also been used to describe neutrino 
propagation through the Earth [221 EH EH]- We will again treat s— and t— channel dominated 
scenarios separately; we saw in the previous Section that this corresponds to scenarios with 
bino- and higgsino-like LSP, respectively. 



3.1 Transport equation for s— channel scattering 

The quantity we wish to compute is the differential flux F$o as function of the LSP energy 
E and matter depth X, where the latter quantity is customarily given as a column depth, 
measured in g/cm 2 or, in natural units, in GeV 3 ; for the Earth, X e [0, 2.398 ■ 10 6 GeV 3 ] |3*U]. 5 
A UHE LSP interacting with matter at rest can only lose energy. The interactions of LSPs 
with energy E will therefore always reduce F^o(E). The size of this effect is determined by the 

5 An LSP exiting the Earth at angle 9 relative to the vertical direction will have transversed a column depth 
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interaction length A^o, which is given by 

X^E)- 1 = N A a^ N (E) , (22) 

where Na = 6.022 x 10 23 g -1 is Avogadro's number, and in the given scenario 0"-°^ — 0"!j ot of 
Eq.fjnj). On the other hand, F^o(E) can be increased by the interactions of LSPs with energy 
E y = E/(l — y) > E losing a fraction y of their energy. These considerations lead to the 
transport equation 

dF^(E,X) F^o(E,X) l /-Km" dv 

^TX- 2 = -i^T + A^(l)i T^ y ^v)F^E v ,X), (23) 

where the kernel is determined by the differential cross section given in Eq.(JHJ): 

K(Ev) = 1 das{Ey) (24) 

and y m ax has been given in Eq.fjHJ). Note that Eq.(J2HJ) assumes collinear kinematics, where 
the produced LSP goes in the same direction as the original one. This is justified for ultra- 
relativistic LSPs, whose scattering angle is 9 < rrig/E in < 1CT 3 for energies of interest. 

An important property of x? — nucleon scattering is that it always produces another Xi m 
the final state. The total flux, 



= / F%(E,X)dE, (25) 



must therefore remain constant, independent of X; here E mSjX is the maximal x? energy, beyond 
which the incident LSP flux vanishes. This is reflected in the fact that integrating the right- 
hand side of Eq. ()2Hj) over the energy gives zero, i.e. d^^o/dX = 0. 6 This allows an important 
check of the procedure used to solve the transport equation. 

The standard method to solve the analogous transport equation for neutrinos is based on 
an iteration [221 El [2H] • It uses an effective interaction length A^o defined by 



F^o(E, X) = i^o(£,0)exp 



X 



A^(E,X) 



(26) 



which in turn leads to the definition of a "Z— factor r 



X 9 o(E) 

W E - X)= l-li(E,X) - (27) 
The transport equation can then be re-written as an integral equation for Z^o. This equation 



can be solved iteratively, starting from the 0— th order ansatz Z: Q = 0. This treatment was first 

Xi 



,(0) 
Xi 

introduced in ref. [22] to describe the propagation of electron and muon neutrinos. In this case 

6 To see this, one re-writes the double integral over E and y into an integral over E y and y and uses the 
definition (122H of the Xi interaction length. 
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the iteration converges fast, so that the first non-trivial solution Zjt} deviates by at most 4% 
from the final answer. The same first-order solution was then used in refs. |2~T] l2~3~] to describe 
v T propagation. We therefore tried to use it for LSP propagation as well. 

Unfortunately we found that the first-order solution badly violates flux conservation in our 
case. A close reading of ref . |22] shows that this is actually not very surprising. To begin with, 
the 0— th order solution evidently ignores regeneration completely. Flux conservation requires 
that the differential flux must increase for some range of energies; Eqs. ()26j) and ()27j) show that 
this is possible only if the Z— factor exceeds unity, which means that the effective interaction 
length A must become negative. This did not happen in ref. [22], since charged current reactions 
of electron or muon neutrinos effectively lead to a loss of this neutrino 7 , thereby reducing the 
total neutrino flux. Flux conservation was therefore of no concern in ref. [22] • Moreover, the 
neutral current reaction, whose regeneration effect was included in this treatment, has da/dy 
quite strongly peaked at low y, see Fig. [7] According to ref. [22], this accelerates the convergence 
of the iteration. However, in our case da s /dy is flat, see Eq.(jHJ). It is therefore not surprising 
that this algorithm does not work very well in our case. 8 

We therefore use a straightforward numerical solution of the transport equation (]2"3~]) . based 
on the first order Taylor expansion: 

F% (E, X + dX) = F% (E, X) + dX ^ + ■ ■ • . (28) 

We found that a more sophisticated algorithm, e.g. the Runge-Kutta method, does not offer 
much of an advantage in terms of accuracy achieved for a fixed CPU time spent. We parame- 
terize F^o(E) for given X as cubic spine. The one-dimensional integral appearing in Eq. (j2rSj) 
can then be evaluated using, e.g., the Simpson method. 

Since the transport equation is a first order differential equation, we have to specify the 
boundary condition F^o(E, 0). We use three different ansatze: 

Spectrum 1:F~ (£,0) = NiE~ 2 F cut (E) ; 



Spectrum 2:F-° (£,0) = N 2 E~ 1 - b F cnt (E) 



I 

(l + lQs'Gev) 



Spectrum 3:F~° ; (£,0) = A^T 1 - ^^F cut (E) , (29) 



where the N{ are arbitrary normalization factors. The first two ansatze would yield an infinite 
total energy stored in LSPs, had we not introduced a cut-off function 

F cut (E) = 1 - (-^-) , (30) 



E 



cut 



7 In the energy range of interest, muons will come to rest before decaying. The produced in this decay 
has such a low energy that it is effectively invisible to neutrino telescopes. 

8 We suspect that the first order solution of this algorithm also leads to sizable errors when applied to v T 
propagation. Since r leptons in the relevant range of energies decay before they interact, the total v T flux should 
be conserved. Moreover, since only a relatively small fraction of the energy of the original v r goes into the v T 
produced in r decay, the effective y distribution, computed as in Ea. (j.S2fl below, should extend to quite large 
values of y. 
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with F cut (E) = for E > E cut . In top-down scenarios, E cnt should be ~ Mx/2. In our 
numerical examples we instead use E cnt = 10 9 GeV. We will see below that LSPs with initial 
energy E > 10 9 GeV will interact very soon, i.e. the flux at E > 10 9 GeV will quickly become 
negligible. Due to the shorter interaction length, choosing E CVLt > 10 9 GeV would greatly 
increase the numerical cost of the solution. On the other hand, this relatively small value of 
E cut -C Mx is justified if neutralinos with E > E cut contribute negligibly to the total LSP flux 
impinging on Earth, which is often the case jTTj. 
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Figure 8: Modification of our three LSP spectra due to propagation through the Earth 
for scenario Dl with bino-like LSP. The solid curve shows the initial spectra (in arbitrary 
units), and the dotted, dot-dashed and dashed curves show the flux after a column depth of 
0.09V max , 0.45X max and V max = 2.4 ■ 10 6 GeV 3 . The top-left, top-right and bottom frames are 
for Spectrum 1, Spectrum 2 and Spectrum 3 of Eqs.([29|). respectively. 



In this Subsection we focus on the first scenario of Tabled We saw in Fig. |2]that this leads 
to the largest cross section, and hence the largest propagation effects. Results for the three 
initial fluxes of Eq. (J29|) are shown in Fig.|HJ The initial spectrum should get reduced significantly 
once the traversed column depth exceeds the interactions length; note that according to Eq. ( 



_tot 



<7~oV = 1 nb corresponds to an interaction length A^o 



0.36 • 10 6 GeV 3 ~ 0.15X n 



Fig. El 
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shows that in scenario Dl this cross section is reached at E{ n ~ 10 7 GeV. We see in the two 
upper frames of Fig. ED that the suppression of the LSP spectrum indeed begins to be noticeable 
around that energy for X = 0.09X max . 

In the lower frame the suppression starts at somewhat higher energy for this small column 
depth, since for this very hard initial spectrum the regeneration effect is maximized. This 
effect creates a bump in the spectrum. The peak of the bump is at an energy such that the 
interaction length is slightly larger than the traversed column depth; since the cross section 
rises with energy, the bump shifts to lower energies with increasing column depth. Conservation 
of the total flux (|29|) implies that this bump is more noticeable for a harder initial spectrum. 
An analogous argument explains why the energy where the spectrum begins to be suppressed 
becomes smaller with increasing column depth. This can also be understood from the simple 
observation that increasing column depth increases the probability for interactions between 
LSPs and matter. 

Fig. [21 shows that at high energies, the total s— channel cross section for bino-like LSP is 
roughly proportional to (l/m^)(Ei n /mg) 0A . According to Tabled squarks are about 2.7 times 
heavier in scenario D3 than in scenario Dl. The onset of the suppression of the initial LSP 
spectrum, and the location of the bump produced by regeneration, will thus occur at about 12 
times higher energy than in scenario Dl. 

As discussed at the beginning of this Subsection, the total LSP flux (|2"5|) must be conserved. 
We checked that our numerical solution satisfies this constraint very well, with maximal devi- 
ation of less than 0.1% even for the hardest spectrum and largest column depth. In contrast, 
for the hardest incident LSP spectrum the prediction of the first-order iterative solution can 
violate flux conservation by a factor of two or more. 

3.2 Transport equation for t— channel scattering 

We consider the transport equation for t— channel scattering only for a scenario with higgsino- 
dominated LSP; Fig. H] showed that t— channel contributions are negligible for a bino-like 
neutralino. As in case of neutrinos [221 EU E3]- there are both charged and neutral current 
contributions. We saw in Fig. El that these predominantly lead to the production of xt an d 
respectively. Here we will ignore the contribution from all other final states; this approximation 
is good at the few % level. 

In case of v T propagation j2H 122] one usually writes two coupled transport equations, for 
v T itself and for the r leptons produced in v T charged current reactions. Correspondingly we 
would need three coupled transport equations in our case, describing the spectra of Xii Xi an d 
X2- However, these heavier produced particles decay well before they lose a significant fraction 
of their energy through interactions. This is true even for r leptons of the relevant energy; 
the xt an d X2 states in the scenarios of interest have much shorter lifetimes [by a factor of 
order (w^ out — m^o) 5 /m^], and significantly shorter decay lengths for a given lifetime (due to 
their larger masses, i.e. smaller 7 factors), t— channel scattering can thus be treated through an 
integration kernel in the transport equation which is computed from a convolution of the x ou t 
production and decay spectra. We treat charged and neutral current contributions separately: 



dF^o(E,X) 



+ 



1 



dy 




(E,y)F^o(E y ,X) 



dX 



1-2/ 
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where the interaction length A x o has been given by Eq.([22]). The structure of the equation is 
very similar to Eq. (j23)) describing s— channel scattering. The integration kernels are most easily 
written as convolutions in the variable z = 1 — y = E ont /E in : 



~nc,cc (w ^ 1 l Zl ' max a a u ,da NC > cc (E y , Zl ) ldT( Zl E y ,z 2 ) 



K^(E,y) = — / dz x dz 2 5(z-z lZ2 ) 



< E ) A 1>min 4 imin ' ■ " dzi T dz 2 

1 f z ^dz 1 da NC ' cc (E y ,z 1 ) idTizxEy^) _ 

"\Z ^minj^^max 



a t (E) J z zi dz x T dz 2 



y=l-z 



y=l—z 

(32) 



The limits 2 1)min max for the outer integration in the first line of Eq.([32|) follow from Eq.()2U|l 
with x — *■ 1, i.e. s — > s = 2E y m^ + m 2 ~ . The limits for the inner integration follow from 
X ou t decay kinematics described below. The functions in the second line of Eq. (j3~2j) ensure 
that y lies within the kinematical limits, with z min = zi iiain z 2:Tnin and z max = z^max-^max- Note 
that both integration kernels in Eq.()31jl are normalized to the total x\~ nucleon scattering cross 
section, which is here approximated by the t— channel contribution given by Eq. ljlUJl . Finally, 
Eq.(|31jl again assumes collinear kinematics, where the LSP produced in x out decay goes into 
the same direction as the original LSP. 

The missing piece in Eq. fj32|) is the differential decay spectrum of the produced Xout- Due to 
the small Xout — Xi m ass difference for higgsino-like LSP, see Table El we only need to consider 
three-body decays, x out — > Xifih, where fi, f 2 are two SM (anti)fermions whose masses we 
neglect. In the x ou t rest frame we then have: 

(TV 1 f^maz 

1 dE* f \M\ 2 , (33) 



dE% 8m Xout (2vr) 3 J E} ^ 

where E% is the energy of one of the two massless (anti)fermions (the energy of the other being 
determined by energy conservation). The integration limits in Eq. (j33|) are given by 

2rriy.E* — m 2 - — m% 

= X \ X = ■ (34) 

2E; 9 - 2m Xout ± 2 A E%-ml„ 



The total Xout decay width T appearing in Eq.(|3*2"|) can be obtained by integrating Eq. (j3l?j) over 
EZ , the lower and upper integration limit being given by m^o and (ra| + m| out )/(2m^ out ), 
respectively. 

In order to boost into the frame where Xont has energy energy Ey. out , we have to know the 
angular distribution of the produced Xi m the Xout rest frame relative to the Xout flight direction. 
Here we assume an isotropic distribution, which is appropriate for an unpolarized Xout, and also 
if Xout decays via a pure vector coupling; we saw in Subsec. 2.2 that the relevant couplings are 
indeed almost pure vector couplings in scenarios with higgsino-like LSP. This yields 

* r^^jn^ , r-«W- (35) 



dE x\ Je: . Xl ^xout p. [E» - m\ 16m Xout (2tt) 



E /,min 
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The limits for the inner integration in Eq. ()35j) have been given in Eq. (j34jl . and the limits for 
the outer integration are: 



E*~ . = <y(E*o - BjE\ -m 2 ) 



171% + m\ 

^o max = (36) 

Xout 



Here, and in Eq. (J3H|) . 7 = 1/ \J\ — 3 2 = Ey, out /m^ out . Notice that E* min reduces to the absolute 
kinematical minimum of m^o for E^o = 7171^0, whereas -E^ax is always determined from the x out 
decay kinematics, independent of E^o. Finally, in the relevant limit 7 > 1 the limits on the 
energy E^o of the LSP produced in Xont decay are 



-£/,~o 



7 



m~ + mi — 3 I mi — m~ 

X" Xout r V Xout Xi 



m 2 ~o 
xi 

7 



E^o 



7 



Xout Xout 



m.Q + mi + 3 I mi — m~ 

xi xout ^ v xout xi 



«' max 2m, t 

Xout 



7 m xout = E xout ■ (37) 



Since the heavier higgsinos Xout have very similar masses as the LSP, see Table El the energy 
loss in the Xout decay process cannot be very large, and may be zero. For example, in scenario 
H2 the maximal energy loss is 6.7% in x° decays produced in neutral current events, and just 
2% in xf decays produced in charged current events. For simplicity we therefore treated these 
decays with purely kinematical decay distributions, i.e. we used a constant matrix element Ai 
in Ens.fPj) and 

Our results for scenario H2 with a higgsino-like LSP are displayed in Fig. El Although 
comparison of Figs. 0] and El shows that this scenario and the scenario Dl with bino-like LSP 
lead to similar total Xi~ nucleon scattering cross sections, comparison of Figs.|Sland|nishows that 
the spectrum of higgsino-dominated LSPs is modified much less by propagation effects. The 
reason is that the t— channel exchange contributions have a far steeper y— distribution, cf. Fig. [3 
The spectrum will be modified significantly only once a given LSP loses a sizable fraction of 
its energy through interactions. The characteristic length for this is better characterized by 
A|o = A^o / (y) than by the interaction length A^o itself. Here (y) is the average relative energy 
loss in an interaction and, for the t— channel, subsequent x out decay; numerically, (y) ~ 0.5 
for scenario Dl, whereas (y) ~ 0.07 for scenario H2. A-o therefore differs by about a factor 
of 10 between scenarios Dl and H2. As a result, even after traversing the whole Earth in 
scenario H2 only the spectrum at E > 10 8 GeV is suppressed, whereas in scenario Dl with 
bino-like LSP the suppression already starts at E ~ 10 6 GeV. The same effect also explains 
the reduced size of the "regeneration bump" in the spectrum of higgsino-like LSPs. Indeed, for 
the steepest input spectrum (top-left frame) this feature is no longer visible; however, Table |U 
shows that flux conservation is satisfied also in this case. 9 We should mention here that due 



9 Due to the steeper y distribution, the first iteration of the method developed in ref.[22] works slightly better 
in this case than for bino-like LSP. However, after traversing the whole Earth we still find that up to 40% of 
the original LSP flux is "lost" if this method is used. 
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to the additional integrations required in Eqs. (|H2|) and (jSSJ), this scenario is numerically much 
more time consuming than the case of a bino-like LSP. 
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Figure 9: Modification of our three LSP spectra due to propagation through the Earth for 
scenario H2 with higgsino-like LSP. Conventions are as in Fig. El 



4 Summary and conclusions 

The detection of ultra relativistic cosmic neutralinos xl would be a crucial test of both "top- 
down" models and supersymmetry. This should in principle be feasible by using the Earth to 
filter out very energetic neutrinos, which would otherwise form an unsurmountable background. 
In order to assess the prospects of future or planned experiments such as EUSO [3T] and OWL 
[32J a realistic treatment of the interactions of neutralinos with ordinary matter is essential. 
In this paper we provide two important ingredients of this treatment: an improved calculation 
of the LSP-nucleon scattering cross section, and a reliable treatment of propagation effects by 
means of transport equations. 

The calculation of the most important contributions to the x° — nucleon scattering cross 
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^/-^max 


Spectrum 1 


^o(X)/^o(0 

Spectrum 2 


) 

Spectrum 3 


0.09 
0.45 
1. 


1.001 
1.001 
1.001 


1.000 
1.001 
1.002 


1.003 
1.013 
1.022 



Table 4: Test of conservation of the total LSP flux (|25|) for the three incident spectra defined 
in Eqs. (j29|) using scenario H2 with higgsino-like LSP; the LSP spectra are integrated between 
10 3 GeV and E cnt = 10 9 GeV. The column depths are in units of the maximal Earth column 
depth X max = 2.4 ■ 10 6 GeV 3 . 



sections are presented in Sec. 2. New features include the calculation of the cross section differ- 
ential in the energy of the outgoing neutralino, which is needed for the quantitative description 
of neutralino propagation through the Earth, and a first calculation of t— channel scattering 
processes, which dominate in case of higgsino-like LSP but are negligible for bino-like LSP. We 
also calculated the contribution from x\ + 9 —> t + 1, see the Appendix. This process turns out 
to be sub-dominant in the scenarios we studied, but might be important in special scenarios 
where the lighter scalar top t\ is much lighter than the other squarks. We found that the total 
cross sections for bino- and higgsino-like LSPs are comparable to each other if first and second 
generation squark masses are around 400 GeV. For lighter (heavier) squarks the cross section 
for bino-like LSPs will be enhanced (suppressed). The cross section for higgsino-like LSP is 
much less dependent on the values of various parameters appearing in the SUSY Lagrangian, 
as long as the LSP is indeed higgsino-like; at high energies it is only a factor ~ 1.5 below the 
cross section for neutrino-nucleon scattering. 

However, the shape of the distribution in the energy loss variable y is quite different in the 
three cases. Under the assumption that squarks are lighter than gluinos, a bino-like LSP has 
a flat y— distribution, with average (y) ~ 0.5. In case of neutrinos, the y distribution peaks at 
y = 0, but falls only gently as y — > 1, with (y) ~ 0.25 |3UJ|2H]. F ma Uy> higgsino-like neutralinos 
have a rapidly falling y— distribution, with (y) ~ 0.07 even after including the energy lost in 
the decay of the heavier neutralinos and charginos that are produced in this case; note that 
t— channel scattering predominantly produces the other higgsino-like states, whose masses are 
however very similar to that of the LSP, so that the energy lost in the decay only amounts to a 
few %. If squarks are lighter than gluinos, most produced squarks decay into a gluino (rather 
than an LSP) plus a quark, with the gluino decaying into at least three particles (including an 
LSP again). This would lead to significantly larger energy loss, i.e. (y) would be even larger. 
This case has yet to be analyzed in detail. Here one could use the formalism we developed to 
treat the decay of the heavier higgsino-like states, although in case of gluinos the decay might 
proceed over several intermediate states |28] . 

The transport equation describing LSP propagation through the Earth are discussed in 
Sec. 3. We first tried the first order iterative solution of ref. [22], but found that it leads to large 
violations of flux conservation in our case; we suspect that this is also true for the rather similar 
case of v T propagation, where the same method was used in [2UI2B]- We therefore solved the 
transport equation by straightforward numerical integration. We found that the y— distribution 
is crucial for understanding the effects of propagation. Indeed, the product of cross section and 
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(y) is a much more reliable "figure of merit" here than the cross section itself. In particular, 
we saw that propagation effects in a scenario with higgsino-like LSP are much milder than in 
a scenario with bino-like LSP, which has similar total cross section a but much larger a{y); in 
the former case the LSP flux is suppressed only for E > 10 8 GeV, whereas in the latter, the 
onset of suppression is at E ~ 10 6 GeV. Of course, the smaller (y) of higgsino-like LSPs also 
implies that the visible energy released by LSP scattering in or near a detector is much smaller 
than for a bino-like LSP of equal energy. The visible spectrum for higgsino-like LSPs emerging 
vertically out of the Earth would nevertheless extend to at least an order of magnitude higher 
energies than that of bino-like LSPs, for similar total cross section. 

Note also that the effective (y) for u T interactions, including the effect of r decays, is about 
0.55, similar to the case of bino-like LSP. Hence propagation effects will deplete the v T spectrum 
much faster than the spectrum of higgsino-like LSPs. We therefore expect that UHE higgsino- 
like LSPs, which have not yet been considered in the literature, could also yield a viable signal 
for visible energies roughly between 10 5 and 10 7 GeV. However, a detailed investigation of 
signals is beyond the scope of this paper. In the analysis of possible signatures for ultra- 
relativistic higgsino-like LSPs, the finite flight path of the heavier higgsino-like states may also 
have to be taken into account for scenarios with mass splitting to the LSP below ~ 2 GeV. 

In this paper we did not discuss the possibility of a wino-like LSP, as in scenarios with 
anomaly mediated supersymmetry breaking J5|. In this case we expect a strongly suppressed 
t— channel neutral current contribution, since a neutral wino does not couple to Z bosons, but 
somewhat enhanced charged current contribution, since the wino is an 577(2) triplet rather than 
a doublet. However, the produced chargino may now lose a significant amount of energy before 
decaying, since the mass splitting to the LSP if typically only a few hundred MeV in this case, 
leading to a lifetime which is much longer than that of the r lepton. The propagation of these 
charginos would then have to be described by a second, coupled, transport equation. A wino- 
like LSP can also have a sizable s— channel contribution to the cross section, if squarks are not 
too heavy. In this case left handed, SU(2) doublet, squarks would be produced predominantly; 
if lighter than gluinos, these squarks would decay into charged or neutral winos, i.e. into the 
LSP or the lighter chargino. 

In summary, we provided a treatment of the interactions of ultra-relativistic higgsino- or 
bino-like neutralinos with matter, calculating both the total cross section and the energy lost 
in these interactions. We also developed a new method to solve the transport equations de- 
scribing the passage of these neutralinos through matter. This will allow improved estimates of 
the prospects of future experiments to observe these neutralinos in potentially realistic super- 
symmetric scenarios. Such an observation would help to solve the forty year old puzzle posed 
by the cosmic rays at ultra-high energies. 
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A Top + stop production 



t,qi + t,q 2 i} 




M, 



M, 



Figure 10: Feynman diagrams contributing to Xi + 9 ~ > t + 1. The arrows indicate the flow of 
baryon number. The diagrams for Xi + 9 ~ > t + t can be obtained by reversing the directions 
of these arrows. 



The Feynman diagrams for top + stop production via x° — gluon fusion are shown in Fig.fTUl 
Here, Pi, kj are four-momenta in the initial and final state, respectively, and qi = P2 — k 2 and 
Q2 — Pi — k 2 are the exchanged four-momenta in the two diagrams. The total squared partonic 
production amplitude is then given by \M a + Mb\ 2 = \M a \ 2 + \Mb\ 2 + M a M* b + M* a Mb- In 
the following we list these three contributions separately. 



\M a \ 2 term: 



spins colors 



\Mb\ 2 term: 



2 (t - ml) 2 



- I 2 + \ G f I 2 ) (fa -qiPi-qi-h- Pi qi-qi + h- q\ qi ■ Pi) 

+ 2G L i Gf*m^m t q\ + 2Gf *Gf 'm^m t q\ 
+ mGf\ 2 + \G«\ 2 )m 2 t p 1 .q 1 
- 2Gf*Gfm^om t k 2 - qi - 2G^Gf*m^m t k 2 -q 1 
-(|Gf| 2 + |Gf| 2 )m^-fc 2 

+ 2G^Gfm^ml + 2Gf G?*m^mf\ . (3c 



spins colors 



u — m~ 



(t + 2k 2 -q 2 + ml^ 



\Gf | 2 + |Gf | 2 ) pvh + 2G L i Gfm^m t \ 



(39) 
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M a M* b + M* a M b terms: 



E (Mow + M * M * 



spins colors 



2 (t — m?)(u — m 2 ~ ) 



\Gf\ 2 + \G^) 



(pi-k 2 q v q 2 +Pi-qik v k 2 -pi-k 1 k 2 -qi + Pi-k 2 k 1 -q 1 



mf mi 
-Pi-to-Y -Pvk 2 — 

+ (<?r<?2 J rk v qi-k 2 -q 2 - k x -k 2 ) (m t m i oGfG i 



(40) 



Here, g s is the strong (QCD) coupling constant, t = q\, u = q 2 , m t is the mass of the top quark 



L Ft 

and mi that of the produced stop squark. The couplings G i ' depend on whether the lighter 
(i = 1) or heavier (i = 2) stop squark is to be produced: 



G\ = -V2g 2 (^~Nt 2 + ~ta,n6 w NZi)cos6 



G 2 



G* 



92m t . 
\/2M w sin p 
92m t 



V2g 2 (\n* 12 + ~ tan 9 W N* U ) sin 6~ t - - N* u cos 6~ t , 

2v/ 2 + a at - a 92m t 
g 2 tan uw-^u sm "t 7= ; -^u cos 5 



3 

2\/2 



g 2 tan 9wNu cos 6> t ~ + 



v^Mh/ sin /? 
g2m t 



where % is the stop mixing angle (i.e. ti = cos 9 ft l + sin %£r). 



7V 14 sin t - , 



(41) 
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Figure 11: The total cross section for Xi + g — * ^> tt (solid curves) in comparison with the s- 
(dotted) and t— channel (dashed curves) contributions. 



The added cross section of Xi9 ~ * tt and Xi9 ~ * tt for a higgsino- and bino-like neutralino 
are displayed in the left and right frames of Fig. ^2 respectively; we used = 1.2 TeV, 
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m t ~ 2 = 1.7 TeV, &i = 0.03 in scenario H2, and = 0.76 TeV, m t ~ 2 = 0.91 TeV, &i = 1.4 in 
scenario D3. We see that in both cases t + t production only make sub-dominant contributions 
to the total cross section, although in scenario H2 it dominates the total squark production 
cross section; the reason is that higgsinos couple proportional to the mass of the quark, as 
shown by the second term of each of Eqs. (j41j) . Note also that the average energy loss in 
t + t production will be even larger than in s— channel reactions involving light quarks. These 
processes can therefore be important for higgsino-like LSP if 1m~ tl < < \M±\, IM2] (which 
does not happen in mSUGRA |15j). or in scenarios with bino-like LSP if t\ is much lighter 
than first or second generation squarks (which can happen even in mSUGRA, if \Aq\ is large). 
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